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a: 

A method is presented for reducing general two-loop self-energies to standard scalar integrals 
in massive gauge theories with special emphasis on the electroweak Standard Model (SM). We 
develop a technique for treating the tensor structure of two-loop integrals appearing in self- 
energy calculations. It is used together with the symmetry properties of the integrals to obtain 
a result in terms of a small number of standard scalar integrals. The results are valid for 
arbitrary values of the invariant momentum p 2 , all particle masses, the space-time dimension 
D and the gauge parameters £j (i = j,Z,W). The algebraic structure of the results clearly 
displays the gauge dependence of the considered quantities and allows to perform very stringent 
checks. We explicitly verify Slavnov- Taylor identities by calculating several thousand Feynman- 
diagrams and adding them up algebraically. As an application of our algorithm we calculate the 
light fermion contributions to the two-loop gauge boson self-energies of the electroweak SM. We 
study their gauge dependence and discuss the occurring standard integrals. 
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1 Introduction 



The e + e~ colliders LEP100 and SLC started a new era of precision measurements 
which allow to test the electroweak Standard Model (SM) on its quantum level. In order 
to match the experimental precision radiative corrections have to be incorporated into 
the theoretical predictions. During the last years many calculations of electroweak virtual 
corrections have been carried out at the one-loop level. For the high precision experiments 
at LEP100, however, first order corrections alone are inadequate. Leading second order 
effects are often taken into account by means of renormalization group methods, but 
rather limited results have been obtained for irreducible virtual two-loop corrections. 

In the study of these contributions the self-energies play a central role. These so-called 
oblique corrections are universal, i.e. process- independent, in contrast to the process- 
specific contributions due to vertices, boxes and bremsstrahlung. A number of au- 
thors [2j recently stressed the importance of the oblique corrections for the analysis 
of precision experiments and in reformulations of electroweak radiative corrections using 
effective lagrangians. In view of the measurements at LEP100, in particular a precise 
calculation of the Z-boson self-energy is of interest. 

In the SM so far no complete calculation of a two-loop self-energy has been carried 
out. This fact is due to the complicated structure and large number of the Feynman 
diagrams contributing at the two-loop level. Results were obtained treating the limiting 
cases of a heavy fermion doublet || , a heavy top quark || and a large Higgs mass || . The 
result given in H was used in || for studying the resummation of effects due to fermion 
doublets with large mass splitting. 

In view of the fact that the mass expected for the top-quark is of the order of the 
heavy gauge boson masses and that almost no restrictions can be imposed on the Higgs 
mass a calculation allowing for general values of the top and Higgs mass, the gauge boson 
masses and the invariant momentum p 2 might be of interest. 

In this paper we present a systematic way for treating all two-loop self-energies. It 
is applicable for gauge bosons, scalars, gauge boson-scalar mixing and fermions. In our 
discussion we focus on the electroweak SM, but the method is valid for any renormalizable 
model. The strategy we have adopted is to reduce the amplitudes algebraically as far as 
possible. We develop a technique for the tensor decomposition of two-loop self-energy 
integrals. This generalizes the results for one-loop integrals worked out by Passarino and 
Veltman J!]]. Using this technique and the symmetry properties of the integrals we obtain 
a result in which the Feynman amplitudes are given in terms of standard scalar two-loop 
integrals. This is in analogy to one-loop calculations where the Feynman amplitudes are 
expressed in terms of the standard integrals A Q ,B ,Co and D defined in ||. We will 
show that the reduction to the specified class of integrals is possible for every two-loop 
self-energy. 

The calculations are performed for arbitrary values of all particle masses, the invariant 
momentum p 2 and the space-time dimension D. We work in a general .R^-gauge specified 
by one gauge parameter £j (i = 7, Z, W) for each vector boson. This allows to study the 
gauge dependence of the calculated quantities. The gauge dependence of basic electroweak 
corrections has recently found considerable interest and was studied at one-loop order by 
a number of authors 0, [], [TU], |TT[. At the two- loop level this issue is of even greater 
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importance since as long as the calculation of a complete process, which necessarily is 
gauge invariant, is out of reach it is crucial to know the gauge dependence of the results. 
We will show that our method of representing the results in terms of two-loop standard 
integrals very clearly displays the gauge dependence of the considered quantities. It can 
be read off directly from the algebraic result. There is no need for using an explicit 
analytical or numerical expression of the standard integrals. 

We show that our algorithm for treating the two-loop self-energies is well suited for 
calculations involving a large number of Feynman diagrams. The results consist of a 
relatively small number of standard integrals which can conveniently be studied for further 
evaluation. In contrast to that a direct evaluation of the tensor integrals would in general 
involve a large number of integrals, each one to be treated separately. The results we 
obtain are transparent and have the benefit that very stringent checks can easily be 
performed on them. In analogy to the investigation of the gauge dependence we can check 
Slavnov- Taylor identities directly at the algebraic level. This test is exact, i.e. free of any 
numerical uncertainty. We implemented our algorithm in the computer-algebra program 
TwoCalc [I2|. In order to check its reliability we explicitly verify Slavnov- Taylor identities 



for the self-energies of the 7Z-system involving several thousand Feynman diagrams. 
Every graph is calculated separately and the results are summed up algebraically. 

As an application we treat the light fermion contributions to the two-loop gauge bo- 
son self-energies in the SM. Considering the results obtained in one-loop order the light 
fermions are expected to yield a significant contribution to the complete result. We give 
the results for the 7, Z and W self-energies in terms of two-loop standard integrals for 
general values of p 2 , the gauge boson masses and the Higgs mass. We study the gauge 
dependence of these amplitudes and classify them in several subsets according to their 
behavior under gauge transformations. All standard integrals appearing in the results 
of the light fermion contributions can be solved analytically leading to an expression in 
terms of polylogarithmic functions. We worked this out explicitly. The results will be 



presented in a related paper [13 



The paper is organized as follows: In sect. 2 we classify the two-loop self-energies 
according to their topologies and define the relevant quantities. Sect. 3 is concerned with 
the symmetry properties of the two-loop integrals. In sect. 4 we develop a technique for 
the tensor decomposition of two-loop integrals and show that it is applicable for every 
integral which arises in the calculation of two-loop self-energies. In sect. 5 we consider 
relations used for minimizing the number of occurring standard integrals. Sects. 6 and 
7 deal with the computer-algebraic realization of the algorithm and its use for verifying 
Slavnov- Taylor identities in the electroweak SM. In sect. 8 we give results for the two-loop 
photon self-energy in QED and for the light fermion contributions to the two-loop gauge 
boson self-energies in the electroweak SM. The properties of these results are discussed. 
In the appendix we give a result for the W self-energy and list several reduction formulae 
needed for the calculations performed in this paper. 
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2 Classification of two-loop self-energies 



The two-loop self-energies can be classified according to the topologies of the corre- 
sponding Feynman diagrams. All self-energy topologies which can occur in renormalizable 
gauge theories are shown in fig. |2.1| . We have listed the one-particle irreducible topolo- 
gies, where as usual diagrams containing tadpole lines are included. Furthermore nine 
reducible topologies exist which correspond to products of one-loop self-energies. In the 
following we will call the first eight topologies in fig. |2.1| "generic" two-loop topologies. 

We focus on the unrenormalized self-energies, a proper renormalization can be done 
after the algebraic calculation has been carried out. In order to make the integrals mathe- 
matically meaningful we use dimensional regularization and work in an arbitrary space- 
time dimension D. 

Inserting fields into the topologies and applying the Feynman rules for propagators and 
vertices leads to the Feynman amplitudes. It is convenient to deal with scalar quantities, 
i.e. to begin with a tensor decomposition. For the gauge boson self-energies it reads 

E£f (p) = (-*,«, + ^)£t'V) - Vj f-^\V 2 ) , (2.1) 

from which the transverse part (p 2 ) and the longitudinal part (p 2 ) can easily be 
extracted: 

£tV) = ^ziirsT + ^f)^ip) ; £l'V) = (p) , (2.2) 

p is the external momentum, D the space-time dimension, a, (3 = 7, Z for the (7Z)-system 
and a = W, /3 = W for the W^-boson. 

We write for the mixing of gauge bosons and unphysical Higgs-fields 

E£ 4 (p)=p M E*V)- ( 2 -3) 

Consequently we have 

E Q >V) = ^£f(p) > ( 2 - 4 ) 
p l * 

where i = x,f for the (7Z)-system and the W-bosoia, respectively. 

The fermion self-energies can be decomposed into a vector, an axial vector, a scalar 
and a pseudoscalar part according to 

£(p) =piY> v {p 2 )+ ^75 T, A {p 2 ) + m£ 5 (p 2 ) + m 75 £ P (p 2 ) , (2.5) 

where m is the mass of the fermion and mixing effects have been suppressed for simplicity 
of notation. We obtain 

and 

Z S (P 2 ) = ^Tr(s(p)) ; S P (p 2 ) = ^Trj^p)) . (2 . 7 ) 
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Figure 2.1: The 20 topologies possible for two-loop self-energies in renormalizable gauge 
theories. 
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Adding the inverse of the zeroth order propagator to the self-energies of all orders 
(including the tadpole amplitudes) leads to the truncated one-particle irreducible two- 
point functions. For the Higgs field being a Lorentz scalar this yields 



r H (p) = i(p 2 - m 2 H ) + iH H (p 2 



(2.8) 



and for the Lorentz tensors we use the decompositions specified above. The corresponding 
propagators are obtained as the inverse of these two-point functions. For the neutral gauge 
bosons we have to consider matrices. The transverse part of the inverse propagator matrix 
reads 



D? 1 



(p 2 ^ 



P 



2 - M 2 +^ z (p 2 )) 



(2.9) 



from which the propagators 



D 1 



A? 



Af z 



(2.10) 



follow by matrix inversion. 

We recall that in the i?^-gauges the lowest order gauge boson propagators can be 
written as 



DIM 





- + - 






ft™ 


k 2 — m 2 




k 2 - m 2 /& 




k 2 — mf 



;2.n) 



where i = 7, Z, W. The parameters £j associated with these fields specify the gauge. They 
can be chosen independently. The 't Hooft-Feynman gauge is realized by setting £j = 1, 
£i —>■ 00 defines the Landau gauge while the limit & — >• corresponds to the unitary 



gauge. The unphysical scalar and ghost propagators are of the form i k 2 — mf/£, 

-1 



The 



fermion propagators can be written as + m) 



At 



m 



In the following we will use the shorthand notation 

d D qi f d D q 2 



((•••)) 



i7T 2 (27T/i) 



D-4 



i7T 2 (27T/i) 



D-4 



(...) 



(2.12) 



where q\ and q 2 are the integration momenta of the loop integrals and \x is an arbitrary 
reference mass. 

The Feynman amplitudes we are concerned with can therefore be written as 



2 



ml 



mi 



(2.13) 



where kj is the momentum of the j-th propagator and rrij its mass, j = 1, . . . ,£. The 
m 2 /£j-terms occurring in the propagators given above are simply treated as mass pa- 
rameters. In (|2.13|) it is understood that the masses carry a small negative imaginary 
part. 

The numerator of the Feynman amplitude in general has a complicated structure 
being a function of the two integration momenta, the external momentum p, the particle 
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Figure 2.2: The topology of the "master integral" 



masses and the gauge parameters It will be a central issue of this paper to show how 
all Feynman amplitudes can be reduced to a form where the numerator consists only of 
quantities which are independent of the integration momenta and can therefore be pulled 
out of the integral. The denominator is still of the form ( [2.13| ). This means that the 
amplitudes can be expressed in terms of a class of scalar two-loop integrals 



Ti 



(P 5 



m 1 ,m 2 , 



, m,, 



« 



k 2 



mi 



k 2 



rrir 



7,2 



mi 



(2.H) 



which we call T-integrals. This is reminiscent of the well known result that every one-loop 
amplitude can be reduced to the basic scalar integrals A ,B , C and D |7|, |8j. It should 
be noted, however, that it is by far not obvious that such a reduction is possible in the 
general two-loop case and we only claim it for self-energies. 

We used the double index notation X^...^ to indicate that the subindices of the T- 
integrals refer to the corresponding momenta k^, ki 2 , . . . ki e . The masses are only explicitly 
written as arguments if confusion is possible. If a propagator has mass zero, we indicate 
this with a prime at the corresponding subindex and drop the zero in the list of arguments, 
e.g. 



T 1/2U (p 2 ;m 2 a ,m 2 b ,m 2 c ) = (( 



1 



i- 2 

"2 



mi 



k 2 



mi 



;.2 

4 



mi 



(2.15) 



The dependence on p 2 will be suppressed in the following. 

The topologies listed in fig. [2.1| can also be used to represent the T-integrals. A line 
in the topology carrying zero momentum contributes only a factor (—l/mf) which is 
irrelevant for the scalar integral. The topologies 6 and 7 in fig. |2.1| therefore represent the 
same type of T-integral. The first two topologies in fig. ^7T| are the most general ones since 



s corresponding to the other topologies can be obtained from these by 
-m 2 in the denominator or equivalently shrinking the corresponding 
ogy to a point. The integral corresponding to the first topology, being 
the only one where five propagators with different momenta occur, is sometimes called 
the "master integral" 



all scalar integra' 
omitting factors 
lines in the topo 



12345 



1 



m\ 



ml 



k 2 



mi 



k 2 



mi 



k 2 



mi 



» 



(2.16) 



where the momenta are labeled as indicated in fig. |2.2j . In the following we do not 
manifestly impose momentum conservation, but keep the overcomplete set of momenta 
ki,...,hs. This will be convenient especially for treating the symmetry properties of 
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the integrals and performing the tensor reductions. It is easy to reexpress the momenta 
k\, . . . , he, by the external momentum p and the integration momenta q\ and q 2 : 

h = qi, k 2 = qi + p, k 3 = q 2 - qx, h = q 2 , h = q 2 + p. (2.17) 



With this convention, the integral corresponding to the second topology in fig. |2.1] , for 
example, can be written as Tn234- We will see below that equivalent representations can 
be obtained by certain permutations of the indices. 



3 Symmetries of the two-loop integrals 



As stated above we start with a tensor decomposition (eqs. Q2.1]) , (|2.3|) and Q2.5|) , 
respectively) in order to obtain scalar quantities. The contraction of Lorentz indices, 
reduction of the Dirac algebra and evaluation of Dirac traces can be worked out like in 
the one-loop case. This results in scalar products of momenta (ki ■ kj), (ki ■ p), p 2 in the 
numerator of the Feynman amplitude. The denominator is unchanged. We now implic- 
itly use momentum conservation and express all scalar products as sums of momentum 
squares, e.g. 



{h-p) = \{kl-k\-p 2 



(3.1) 



Subsequently all kf appearing both in the numerator and the denominator are canceled 
via 

kf = (kf-m 2 )+m 2 . (3.2) 

The application of this procedure directly leads to T-integrals, if all kf appearing in the 
numerator can be canceled. If not, we obtain another type of integral which contains 
squares of momenta in the numerator not occurring in the denominator: 



yi... 



(( k i ■ ■ ■ Ajki...)) , i,j,k,l - 
Here we used the shorthand notation 

^ikl = T TT : 

K -m 2 k 



1,...,5, 



i^j,k,l . 



k] - m 2 



kf - mf 



(3.3) 



(3.4) 



We now have obtained a representation of the amplitudes in terms of scalar T- and 
Y-integrals. These integrals, however, are not independent of each other. The number of 
occurring integrals can considerably be reduced by taking into account their symmetries 
with respect to permutation of the kf or, equivalently, of the corresponding indices. 

All T- and ^-integrals are invariant under the permutations 



(12)(45), (14)(25), (15)(24). (3.5) 

Here, as always, k±, . . . , k 5 are labeled according to (|2.17j) . Application of the first per- 
mutation yields for example 

.2 _2\ _ rp i 2.^2 _2 _2 



^"1234 (p 2 \rir,. II) 



^1235(P"; 
(( 



k 2 — m 2 



k 2 -m 2 a 



k 2 



rat 



k 2 



rat 



)) • (3-6) 
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The validity of the symmetries listed in (|3.5|) can most easily be seen for the "master 
integral" (fig. |2.2|) or can be checked using ( |2.17[) and the invariance properties of the 



integrals with respect to changes of the integration momenta. 

Additional symmetry relations hold if an index does not occur in the integral, e.g. 

(23) if 1 is absent. (3.7) 



The symmetry relations Q3.5Q and (3/7) are used to map every integral onto a standard 
representative, i.e. all integrals which are related by symmetry transformations are brought 
to the same form. 



4 Tensor reduction for two-loop integrals 

As claimed above, all self-energy amplitudes can be expressed by an independent set 
of T-integrals. This means that all ^-integrals can further be reduced. To achieve this, we 
rewrite the ^-integrals in terms of tensor integrals and perform a tensor decomposition. 

As a simple example we consider the integral 

^2345 = «^A 23 45» ■ (4.1) 

Insertion of 

k\ = (hi - ml) + {ml + p 2 ) - 2(p • k 2 ) (4.2) 

yields 

^2345 = ^345 + (™2 + p 2 )T 23 45 ~ 2p M ( (A# A 2345 ) ) ■ (4.3) 

Now one has to perform a tensor decomposition for the integral 

S 2 2 ^ = ((^A 234 5» • (4.4) 

It is obvious that the Passarino-Veltman procedure J7| used for one-loop integrals, i.e. the 
ansatz 

= P"S(P 2 ) , (4.5) 

does not lead to simpler integrals in this case. To determine the scalar quantity S(p 2 ) one 
has to contract with p^, but it is not possible to cancel the resulting scalar product (p-k 2 ), 
since it is not expressible as a sum of squared momenta occurring in the denominator of 
the integral. This is due to the fact that the momenta p and k 2 belong to a four-vertex in 
the topology A2345 (see fig. |4.1] ). This is a typical feature of two- loop topologies, whereas 
in one-loop order only three-vertices occur in the loop integral and the scalar products 
can always be canceled. 

In our approach we not only work with the tensor structure of the integral with respect 
to the external momentum p but also perform decompositions with respect to a subloop. 
In general one can write for a subloop s M which has the structure of a first-rank tensor 

s* 1 = rfs 1 +p%s 2 + ... (4.6) 
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where si, S2, ■ ■ ■ are Lorentz-scalars and pi,P2, ■ ■ • are the independent external momenta 
of the subloop. 

At first sight it is not obvious that the scalar products resulting from this decompo- 
sition can be canceled. Moreover, as is well known from one-loop calculations, contract- 
ing ( |4.6| ) with pi l j,,p2 ( i, ■ ■ ■ and solving for s\,S2, ■ ■ ■ leads to factors which are determinants 
of the momenta p^. These, however, depend on the integration momentum of the second 
loop and cannot be pulled out of the tensor integral like in the one-loop case. 

In order to show how this technique can successfully be applied to all two-loop self- 
energies we first examine, which denominators of tensor integrals can occur after the 
steps performed in the previous section. They can be read off from the topologies listed 
in fig. [2.1| . Topology 1 clearly is not possible, since at least one of its five different 
propagators can always be canceled. The topologies 9 through 20 are just products of 
one-loop contributions. For these the Passarino-Veltman technique can be applied in a 
straightforward way. The tensor decomposition for the topologies 6 through 8, having no 
dependence on the external momentum p, is also trivial. 

We therefore have to consider the topologies 2-5. As shown in fig. |47T] they can 
be written as A23455, A234 and A2345, respectively. Note that the scalar integrals corre- 
sponding to topologies 4 and 5 are equivalent due to the symmetry relations ( |3.5| ). The 
topologies we have to deal with all contain at least one subloop which is a self-energy 
insertion, in our terminology this is A23. We do the decomposition for this subloop. In 
the simplest case of a first-rank tensor we have 

(A£A 23 ) = A£ S (A; 5 2 ) (4.7) 

which is just the example considered in ( f4.4| ). k§ is the external momentum of the subloop 
A23. The Lorentz-scalar s(k£) is obtained from 

s(kl) = i<(*5 " fc 2 )A 23 ) = ((h ■ fc 2 )A 235 ,> ■ (4.8) 
% 

Due to the self-energy structure of A23 the determinant 1/k 2 has the form of a massless 
propagator. Insertion of fl4.8|) into Q4.71) and (|4.4j) therefore leads to integrals which still 
belong to the class of T- and ^-integrals. We get 

PA%> = «(P • h){h ■ fc 2 )A 2 3455'» • (4-9) 

Since the momenta k 5 and fc 2 belong to a three vertex in the topology A2345 their scalar 
product is expressible as a sum of squared momenta which can be canceled: 

(k 5 -k2) = \(k 2 2 -k 2 3 + kl). (4.10) 

The same holds for {p ■ fc 5 ) 

(p-h) = l(kl-kl+p 2 ), (4.11) 
if one notes that one of the fc|-terms can be canceled against the massless propagator A5/. 
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*23455 



*234 



*2345 



Figure 4.1: The topologies representing the denominators of tensor integrals for which 
the Passarino-Veltman technique is not directly applicable 

We can therefore express Y 2 \ 45 solely m terms of T-integrals. The result reads 



^2345 = A (ml)B Q (p 2 ;ml,m 2 5 ) - \A Q {m 2 2 ) - A (ml)] \^A (m 



m% 



-B (p 2 ;ml,m 



+ T 2U - T 235 - (m 2 2 + m 2 3 + ml - ml + p 2 )T 2345 



-(ml - mf)T 2355 i + (mj - p 2 ) 
-(m 2 2 - ml)(ml-p 2 )T 23i55 > \ 



T; 



2455' 



T- 



3455' 



(4.12) 



The one-loop scalar integrals Aq and Bq appear since some T-integrals are expressible as 
products of one-loop integrals, e.g. 



where 



T 345 = ((A345)) = (A3XA45) = A (ml)B (p 2 ;mlml) 

1 



A (m 2 ) = ( 



q 2 — m 2 



) 



and we use 



B (p 2 ;ml,ml) = ( 



(...) 



q 2 — m\ 



(q + p) 2 — m\ 



) 



d u q 



(4.13) 

(4.14) 
(4.15) 

(4.16) 



i7r 2 (27r/x) D - 4V ' ' 

in analogy to ( p. 12 ). Note that our definitions slightly differ from those used in ||. 

The integrals containing a massive and a massless propagator with the same momen- 
tum can further be simplified by partial fractioning. 

The topology A23455 is treated in analogy to A2345. The corresponding integrals are 
easily obtained by taking the derivative with respect to m\, e.g. 



Y 1 

1 23455 



d 



-Y 1 



(4.17) 



The decomposition for the tensors of higher rank can be worked out in a similar manner 
as described above. In order to treat ^2345 we nave to decompose the integral 



«(p • k 2 )(p ■ k 2 )A 2U5 )) = PlM p u (A 45 (k%k u 2 A 23 )) . 



(4.18) 
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We write 



(A^Ajh) = sqo^V + (4-19) 

with the Lorentz-scalars s 00 (k 2 ) and Su{k 2 ). Contraction with g^ v and k 5fJj k 5u , respec- 
tively, leads to two equations for the scalars s o and Su: 

s 00 D + su = (A| A 23 ) 

s 00 + s u = -z((k 5 • k 2 )(h ■ k 2 )A 23 ) . (4.20) 

Again, the factors l/^f appearing in (|4.19|) and Q4.20|) can be written as propagators 



with mass zero. Inspection of ( [4.20D and ( (4.18[) shows that the decomposition for the 



integral (|4.1S|) leads to T-integrals and one V-integral Y" 2 | 4 which corresponds to a tensor 
of lower rank. The tensor reduction of this integral will be worked out below, where we 
discuss the integrals with the topology A234. The complete reduction formula for is 
rather lengthy. We give it in the appendix. 

The generalization to tensors of higher than second rank should now be obvious. The 
decomposition of the tensor integrals always leads directly to T-integrals or to tensor 
integrals of lower rank. 

So far we showed that the tensor decomposition works for all tensor integrals with the 
topologies A2345 and A23455. The third topology we have to consider is A234. This case 
is slightly more complicated due to the presence of the two four-vertices in the topology. 

The F-integrals which can occur are Y 2U , F 2 | 4 , ^234 > ^234; ^234; The integrals Y 2U 

and Y" 2 5 34 are not independent of each other but are related through a symmetry relation 
which amounts to a permutation of the masses m 2 and m 4 . The same holds, of course, 
for Y 2 \\ and Kf 4 . 

We begin with the integral F 2 34 an d apply exactly the same steps as for F 2 3 45 . This 
gives 

YL = T34 + (m 2 2 +p 2 )T 2U - 2p^(A^A 23 )) 

= T 34 + K +p 2 )T 23 4 - 2«(p • k 5 )(k 5 ■ k 2 )A 2U5 ,)) , (4.21) 



where ( ^4.7| ) and (O) were used. The last integral looks very similar to (4^) but here 



it is not possible to cancel the k\ term resulting from the insertion of ( |4.10|) and ( [4.11|) 
Instead we obtain 

Y 2 \ 4 (ml ml ml) = ~Y 2 5 u (ml m% m\) + f(T) , (4.22) 

where f(T) summarizes terms in which only T-integrals occur. In order to find an expres- 
sion for I234 in terms of T-integrals only, we use the last symmetry relation listed in ( |3.5| ) 
and obtain after appropriate relabeling of the masses 

^34(^2, ml ml) = - l -Y 2 \,{ml ml m 2 A ) + g(T) , (4.23) 



where g(T) again depends on T-integrals only. Insertion of this equation into ( f4.22|) gives 
the result 

^234 = o\ h ( A o, Bo) + (m 2 2 + m 2 3 + ml+ p 2 )T 2U + 2{m\ - m\)T 2W + (m\ - m\)T V3A 
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+2(mj - ml)(ml - p 2 )T 23Ab , + {m% - p 2 ){m 2 3 - m\)T V23A 



(4.24) 



where h(A , B ) represents a function containing only one-loop integrals. 

Another complication can be seen by inspecting ( |4.24j ). Due to the absence of the 
momentum k 5 the integral Y 2 \ 4 is symmetric with respect to the permutation of the 
momenta k 3 and k A or, equivalently, of the masses m 3 and m 4 . However, the tensor 
decomposition artificially introduces the momentum k 5 into the terms on the right hand 



side of ( 4.24 ) and the symmetry seems to be lost. It can be made manifest again by using 
the relation ( |5.3| ) stated below. The result for Y 2U finally reads 



Y 1 

2 234 



-I A (m 2 2 )A (ml) + A (m 2 2 )A (ml) + A (m 2 3 )A (ml) + (m\ + m 2 3 + ml + p 2 )T 23A 



-{m\-p 2 ) 



A (m 2 2 ) - A (m 2 3 ) B (p 2 -mlO) - (m 2 3 - p 2 ) A (m 2 2 ) - A (m 2 



xB (p 2 ; m 2 3 , 0) + (mj - ml)T 23h > + [m\ - m 2 A )T 23 y{m\, m 2 A ) 

, 2 2w 2 2x , 2 2w 2 2, , 2 2 2^ " 

+ (m 2 - m 3 )(m 4 - p )T 23A y + (m 3 - p )(m 2 - m A )T 23 ^{m 2l m 4 , m 3 ) 



(4.25) 



which displays the considered symmetry. 

The integrals involving higher tensors are treated in exactly the same way. For the 
second-rank tensor we use ( |4.19| ). The reduction of Y 23 \ leads to a formula involving 1^34- 
For this an expression analogous to ( |4.23|) is obtained. The symmetry properties of the 
integral become manifest after using a relation involving T-integrals. The expressions for 
1^234 and are listed in the appendix. 

With the technique for the tensor decomposition of two-loop tensor integrals described 
above all F-integrals can be reduced to scalar integrals of simpler structure where no 
integration momentum appears in the numerator, i.e. to T-integrals. The class of T- 
integrals therefore suffices to express all possible two-loop self-energies. 

This is true in every renormalizable gauge. For practical purposes we briefly examine 
the situation in the 't Hooft-Feynman gauge. As explained above the most complicated 
Y-integrals have denominators of the form shown in fig. |4.1| . The possible numerators can 
be obtained by power counting. For = 1 the gauge boson propagators ( p.ll|) do not 
contribute powers of momenta to the numerator of the Feynman integrals. Counting the 
powers arising from the other Feynman rules and from the decompositions ( p.2|) , 



( |2.6p and ( |2.7p reveals that the relevant integrals are Y 2Ml Y 2 \ Ar)1 Y 23 \ b1 Y 2 \ ibb and Y23455 • 
Explicit formulae for these integrals are given in this section and in the appendix where 
we also list some reduction formulae needed for general values of 

By using the methods outlined above the task of evaluating each Feynman amplitude 
corresponding to a two-loop self-energy is reduced to the calculation of a relatively small 
number of scalar integrals possessing a simpler structure. This is very desirable since a 
direct evaluation of the tensor integrals appearing in the Feynman amplitude by means of 
Feynman parameters in general yields a large number of different integrals which all have 
to be treated separately. In contrast to that the T-integrals are well suited for studying 
their analytical properties as was already done by a number of authors. 

In [14|, [15], |I6|| , for example, results in terms of polylogarithmic functions were 
obtained for special cases of masses and momenta. In more general cases this class of 
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functions is not sufficient. As was shown in | I7| , the general massive two- loop integrals 
cannot be expressed in terms of Nielsen polylogarithms with arguments being algebraic 
functions of the external variables. For the "master integral" (fig. |2.2| ) an integral repre- 
sentation suitable for numerical evaluation was derived in KIT 



5 Relations between the scalar integrals 

The results obtained for the two-loop self-energies via the steps described in the pre- 
ceding sections consist of a Lorentz tensor specified in ( |2.1| ), ( |2.3| ) and ( [2.5| ), respectively, 
and a scalar part being a sum of T-integrals multiplied by rational functions in the squared 
momentum p 2 , the space-time dimension D, the particle masses rrij and the gauge param- 
eters £j. The algebraic structure of the result is convenient for studying the dependence on 
these variables, either by inserting the exact results for the integrals, if these are known, 
or suitable approximations, e.g. in the asymptotic or threshold regime of p 2 or in the limit 
D -> 4. 

The T-integrals occurring in the result, which we denote as T , T 2 , . . . , T n for the mo- 
ment, are in general not algebraically independent of each other, i.e. there exist relations 

c 1 T 1 + c 2 T 2 + ... + c n T n = 0. (5.1) 

The coefficients c 1 , . . . ,c n are polynomials in p 2 ,D,rrij and the gauge parameters. The 
results can be made very transparent by using these relations to eliminate as many inte- 
grals as possible. For example, if one considers a gauge invariant set of amplitudes this 
property directly manifests itself in the algebraic result provided that the set of integrals 
used for expressing the result is sufficiently small. In this case all terms depending on 
the gauge parameters £j disappear from the result. This means that the prefactors of 
all gauge dependent basic integrals algebraically add up to zero while all ^-dependent 
terms multiplying a gauge independent basic integral exactly cancel each other. In this 
way the gauge invariance is seen to hold exactly, i.e. by purely algebraical means. No 
explicit analytical or numerical expressions of the standard integrals are needed for this 
consideration. 

In the same way as the gauge invariance Slavnov- Taylor identities or any other relation 
involving two-loop self-energies can be checked exactly. This is very useful for performing 
consistency checks on the result. 

The question arises whether a sufficiently small basis of integrals can actually be 
found. We argue that after invoking the symmetry properties of the two-loop integrals 
and eliminating all F-integrals by means of tensor decompositions only few additional 
relations between the scalar integrals are needed. In sect. [7] we will verify Slavnov- Taylor 
identities valid for the self-energies of the photon and the Z-boson by adding up the results 
of several thousand Feynman amplitudes. For this application only the relations for the 
integrals ^^/(m 2 ,, m\, m 2 ) and 12345/5/ (m|, m|, m 2 ) described below are used. 

The first type of relations we want to consider involves integrals in which at least 
one propagator is massless. They are obtained indirectly by using symmetry arguments, 
momentum conservation and properties of the tensor decomposition. In this way we 
derive a relation between the integral T^s^m 2 ,, m 2 , m 2 ) and the integrals of the same 
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type where the masses are permuted. It is a consequence of momentum conservation at 
the four- vertices of the topology A234. Starting from 



P ■ (fa + h - k A - p) 



k'2 



rrir 



k 2 



mi 



k\ 



m 



» = 



(5.2) 



which obviously follows from ( |2.17| ), using the symmetries of the resulting integrals with 
respect to permutations of k 2 , k 3 and k 4 and performing tensor decompositions for these 
leads to the relation 

{(m 2 2 - ml)(ml - p 2 )T 23 45'(rnl,ml,ml) + cycl.} 
+ [(m\ - ml)T 235 >(m 2 2} m 2 3 ) + cycl.} 

-{ [A (mj) - A (ml)\ (ml - p 2 )B (p 2 ; ml 0) + cycl.} = . (5.3) 

Here "cycl." denotes cyclic permutation of the masses. With (|5.3|) one permutation of 
the T2345'-integral can always be eliminated. This relation was used in the last section to 
obtain ( f4.25| ). If two masses are equal, ( |5.3| ) becomes trivial. 

Two similar relations hold for the permutations of the integral T 2 u5'5'- Using these, 
it is possible to obtain one standard permutation for every integral of this type. In 
contrast to ( |5.3|) this relation remains nontrivial if two of the three masses are equal. It 
gives an expression for the integral T 234 w(m 2 , M 2 , M 2 ) in terms of integrals with fewer 
propagators only. Therefore it can completely be replaced by integrals of a simpler type. 
We list this relation in the appendix. It takes a particularly simple form in the special 
case where M 2 = 0: 



2~23'4'5'5'( m 



Dm 2 p 2 



A (m 2 ) \(D - 4)B (p 2 ; 0, 0) + Dp 2 B' (p 2 ; 0, 0) 



-(D - 2) 2 T 2 ^(m 2 ) + (3D - 8)T 23/i ,(m 2 ) 
-(D-4)(m 2 +p 2 )T 23 , 4 , 5 ,(m 2 )) . 



(5.4) 



The one-loop integral B' is defined as 

B' (p 2 ; m ii m i) 







d(m 



2 Mp 2 ;mlm 2 2 ) 



(5.5) 



Relations similar to those valid for T2345/ and T2345/5' can be obtained for integrals with 
more massless propagators, e.g. T2345/5/5/ and T 23 4 5 i 5 > 5 r 5 >. 

A second class of relations can be derived via the well known method of integration 
by parts. It yields relations for integrals where at least one propagator appears with a 
power higher than one. A particularly interesting example is the formula for the integral 
corresponding to topology 6 in fig. [2~l] , i.e. Tii 34 (m^, mf, mj, m 2 ). Starting from 



131 



«Al34^» = -«A^A 134 » 



l dk1 



(5.6) 
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and using the symmetry of the integral T 134 with respect to permutation of the masses 
m\ and m| leads to 



(3 - D){m\ -ml- m\)T lu - (1 - D/2) A {ml)(A (ml) 



^1134 ~~ 77 2 2 2^ 

A{mf , m 3 , m|J 



where 



\ / 2 2 2\ 4 , 4, 4 o/22i 22, 2 2\ 

A(m 1 , m 3 , m 4 J = m x + m 3 + m 4 — /(rr^mg + m 1 m 4 + m 3 m 4 ) . 



(5.7) 



(5i 



The integral Tn 34 can therefore be expressed through Ti 34 and products of one-loop inte- 
grals. 

Finally it should be noted that some integrals vanish trivially in the framework of 
dimensional regularization, e.g. 



11 



12 



112 







(5.9) 



6 Computer-algebraic realization 

The procedure we described for the algebraic calculation of two-loop self-energies is 
algorithmic and we implemented it into a computer-algebra program which we called 
TwoCalc fL2| . It is written in Mathematica and linked to the packages FeynArts |1| 
and FeynCalc [|2L)[| . FeynArts creates the Feynman amplitudes and draws the Feynman 
graphs. FeynCalc is used here for the contraction of Lorentz indices and the evaluation 
of Dirac-traces. In the calculations performed in this paper an anticommuting 75 in D 
dimensions was used. 

The three packages run fully automatically. This high degree of automation is very 
convenient for performing calculations which involve a large number of diagrams. 



7 Slavnov- Taylor identities for two-loop self-energies 

In order to demonstrate the abilities of our algorithm for doing large calculations we 
explicitly verify some Slavnov- Taylor identities valid for two-loop self-energies. In this 
section we use the 't Hooft-Feynman gauge. 

In ^T| exact relations for the gauge boson propagators A^, the gauge boson Higgs 



mixing propagators A^ 1 and the unphysical Higgs propagators A 13 are listed 

pVAJzCp) = -1 

p^Ap(p) -iM zP »AJx(p) =0 1 ' 

p ^Af u w (p) + 2M w p»Ajv(p) + M^A^(p) = -% , 
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where \ and (p are the neutral and charged unphysical Higgs fields, respectively. Writing 
the propagators as the inverse of the truncated one-particle irreducible two-point functions 
as specified in sect. ^]and expanding up to second order leads to the second-order relations 



P 2 
P 2 
P 2 
P 2 



£77.(2) _ ( S7X ,(1) 



2 







p 2 E 7x,(l) S Zx,(l) + ^M Z S 7X '( 1 )S»'( 1 ) =0 



The superscripts ^ indicate the two-loop and one-loop self-energies, respectively. Similar 
relations can be obtained for the one-loop self-energies [|2T|. As always, it is understood 
that the tadpole contributions are part of the self-energies (see fig. |2.1|) . 

The first relation in ( |7.2| ) indicates that unlike in pure QED the longitudinal part of 
the two-loop photon self-energy does not vanish. We used TwoCalc to calculate explicitly 
every Feynman amplitude contributing to the left hand side of this relation. We did this 
in the full electroweak SM with the only restriction that the fermions are limited to one 
doublet, which we chose to be the u- and <i-quark. The calculations for all other fermions 
can be performed analogously. 

We had to consider 1024 graphs for the two-loop photon self-energy, 162 of which 
contain the u- and ci-quarks, the remaining ones are purely bosonic. The most involved 
diagrams are of course those stemming from the generic two-loop topologies. These are 
418 diagrams belonging to the first 6 topologies in fig. |2.1| . The remaining ones lead 
to integrals expressible solely as products of one-loop integrals. Altogether 13 of the 20 
possible topologies are realized for the photon self-energy. 

The one-loop (7x)-mixing self-energy vanishes for fermionic graphs, while 4 purely 
bosonic graphs contribute. Summing up the results of all amplitudes gives exactly zero, 
i.e. the prefactor of every standard integral vanishes separately. 

We also explicitly verified the second relation stated in ( |7.2|) . Here, of course, a lot 
more Feynman diagrams have to be calculated. Again we restricted the fermions to u- 
and ci-quark. This gives 1552 diagrams for Ti ZZ '^ 2 \ 1155 for £ Zx >( 2 ) and 1142 for £ xx '( 2 \ 
From these, 1787 diagrams arise from the 8 generic two-loop topologies. 

For the self-energies of the Z-boson and the unphysical Higgs field x au 20 topologies 
contribute which are maximally possible for two-loop self-energies. The self-energies in 
one-loop order give rise to 65 diagrams. Adding those to the 3849 two-loop graphs by 
using the coefficients specified in (|7.2j ) gives exactly zero. 



8 Results for the electroweak Standard Model and discussion 

As an application of the techniques described above we treat the gauge boson self- 
energies in the electroweak SM. We begin with the special case of the photon self-energy 
in pure QED. The three contributing diagrams are shown in fig. |8.1| . 
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Figure 8.1: The diagrams contributing to the two-loop photon self-energy in QED. 



We decompose the photon self-energy according to (|2.1|) . The longitudinal part van- 
ishes due to the well-known Ward-identity. We checked this by explicit calculation. The 
result for the transverse part reads 



£77.(2) ,,;), 



J T,QED 



(P 2 



128tt 4 (D - 1) v 
xA (m 2 e )B (p 2 -0,m 2 e ) 



(2 - Df/ml\Al(ml) - [(2 - D)\m\ - p 2 )/m 

"(2 - D){Am 2 e + 2p 2 - Dp 2 )/ml]A (m 2 e ) 
xB (p 2 ;m 2 e ,m 2 e ) - (8m 2 e - Up 2 + 9Dp 2 - D 2 p 2 )B 2 (p 2 ;m 2 e ,m 2 e ) 
-(2 - D) [2(4m 2 e - 2p 2 + Dp 2 )A (m 2 e )B' (p 2 ; m 2 e , m\) + (2 - D)T m 

-2(1 - D)T 134 , - (6 - D)T 2M , + (2 - D)(m 2 e - p 2 )T V2M , 
+8{4m 2 e - 2p 2 + Dp 2 ) [T 12U , + m 2 e T ll2 - M 

+2{2m 2 e - p 2 )(Am 2 e - 2p 2 + Dp 2 )T 123 , 45 \ . (8.1) 



As above a prime at the subindices of the T-integrals denotes that the corresponding 
propagator has mass zero. All other propagators carry the electron mass m e . 

The calculation was carried out in an arbitrary i?^-gauge. The absence of the gauge 
parameter £ 7 indicates that the result is gauge invariant as it has to be. The two- loop 
integrals appearing in ( |8.1|) are expressible in terms of polylogarithmic functions, i.e. log- 
arithms, di- and trilogarithms. They were studied for example in [15, 16]. The one- loop 



integrals arising from the calculation of irreducible two-loop contributions are in general 
needed up to 0{D — 4). Results for these can be found in |13|, |22| . 

Calculations of the two-loop photon self-energy in QED were performed in ||23|| . Similar 
calculations focusing on the (D(aa s ) corrections due to the gluon exchange in a quark loop 
of gauge boson self-energies were done in [J24|J . 

In passing we note that in contrast to the photon self-energy the two-loop electron 
self-energy in QED is not gauge invariant. The explicit algebraic result will be presented 
elsewhere. It involves the integrals T 23 4,Ti/ 2 34 and Ti/2345' which cannot be reduced to 



polylogarithmic functions. The problem of evaluating Ty 234: 5' was addressed in [[15 



In the electroweak SM results for two-loop gauge boson self-energies restricted to the 
region p 2 = and asymptotic values of the top and Higgs mass, respectively, were worked 
out in [[], [|, |I| . A calculation allowing for general values of the invariant momentum p 2 
and the masses of the gauge bosons, the top quark and the Higgs field has to deal with an 
enormous number of Feynman diagrams having in general a very complicated structure. 
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With the methods outlined above we can express these in terms of a small number of 
standard integrals, the result being valid for general values of all parameters involved. 

This was explicitly carried out in sect. [7] for all graphs contributing to the longitudinal 
parts of the photon and Z-boson self-energies. The results for the complete transverse 
parts are of course rather lengthy and we do not present them here. Instead we focus on 
the light fermion contributions, a class of diagrams particularly well suited for separate 
treatment which is expected to yield a substantial contribution to the complete result. 

With light fermion contributions we mean all diagrams containing any fermion other 
than the top-quark. These fermion masses are negligible in comparison to the boson 
masses (m^,mv[/,m^) and the invariant momentum p 2 , if p 2 is sufficiently large. All T- 
integrals resulting from these graphs therefore contain one massless fermion loop. This 
feature allows an analytic evaluation resulting in polylogarithmic functions. We have 
explicitly calculated all integrals needed for the two-loop gauge boson self-energies, i.e. all 
integrals appearing in the algebraic results stated below. The results will be presented in 
a related paper fll3| . 

In order to study the structure of the results we find it sufficient to restrict the fermions 
to one doublet in the same way as in sect. [7]. This yields 27 diagrams contributing to 
the two-loop photon self-energy, 32 graphs for the Z self-energy and 45 diagrams for the 
self-energy of the W . 

As shown in fig. |8.2| it is convenient to subdivide the photon self-energy into several sets 
of graphs. The graphs depicted in fig. |8.2j a are of QED-type having a photon exchanged in 
a quark loop. As indicated in the picture the same set of graphs also exists for the u-quark. 



The results are in exact analogy to (|8.1|) and we do not write them out explicitly. As 
a consequence of ( |8.1| ) the QED-type graphs form gauge invariant subsets of the photon 
self-energy. 

In the graphs shown in fig. |8.2[ b a Z-boson is exchanged instead of the 7. We give the 
result neglecting the mass of the (i-quark: 

E^jfV) = "'I 2 ~J ){9 ~ US2W t 8 ^ {2(4 - QD + D 2 )A (m 2 z )B (p 2 ; 0, 0) 
T ' z ^ ' (4tt) 4 (1 - D)D\mc 2 w s 2 w \ y J UV ZJ 0Ky ' ' J 

-D(2m 2 z + 7p 2 - Dp 2 )B 2 (p 2 ; 0, 0) + 2(4 - QD + D 2 )T n , 4 ,(m 2 z ) 

+2(8 - AD + D 2 )T 23 , 4 ,{m 2 z ) - 2(4 - 8D + L> 2 )(m| +p 2 )lW 34 ,(m|) 

-D(2m% + 8m 2 z p 2 - Dm 2 z p 2 + 2p A )T VTU , 5 ,{m 2 z ) \ . (8.2) 
The subscript Z refers to the Z-boson exchange. We used the abbreviations 

c w = 2" ) s w = 1 1~ • (8-3) 

m z m z z 

It can be seen from ( |8.2| ) that for non-exceptional values of p 2 no mass singularities are 
induced by setting the fermion mass to zero (see |25| ) . 

There is no gauge dependence left in ( ft.2p . Therefore all contributions from "neutral 
currents", i.e. 7 and Z exchange, are gauge invariant. 
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d <-> u 



(a) 



d ^ u 



(b) 



(c) 

Figure 8.2: The light fermion contributions to the two-loop photon self-energy 
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The "charged currents" mediated by W ± and ^ give rise to the 15 diagrams shown 
in fig. |8.2| c. The result for the transverse part reads 

e 4 (2 - D) 



-<77,(2) (^2 
J T,W± 



(P 2 ) 



f(p ,m w ,D) +9g(£ w ,p } m w ,D) 



(4vr) 4 (l- J D) 2 12m^ L 
+F(p 2 , m 2 w , D) + 9m 2 w {2m 2 w - p 2 )G^ w , p 2 , m 2 w , D) 



■4) 



where in f(p 2 , rriyy, D) and g(£w> P 2 ■, m \yi D) we collected the contributions involving only 
one-loop integrals 



r/40 

f(p 2 , m 2 w , D) = 2m 2 w ( — - 91 + 70D - 10D 



m. 



18(2 -D)p 2 A (m 2 w )B (p 2 ;0,0) 



+8(1 - D)™ 4 , [2m^ + (7 - D)p 2 ] Blip 2 ; 0, 0) - 9(2m 2 w - p 2 ) 



x 



4(1 - D)m 4 w + 4(3 - 2D)m 2 w p 2 - p 4 5 (p 2 ; 0, 0)B (p 2 ; m 2 w , m 2 w 



15) 



9{iw,P 2 ,m 2 w ,D) 



2m 



w , 



(m 



w 



p 2 )A (m 2 w )B (p 2 ;0 } 0) + 2m 2 w 



^ I 2 2\ 
7- (%-P ) 



m 



+2(2 - D)p 2 U o (my£w05 o (p 2 ;O,0) - 2 (™w - P 2 ) \ m w(j- ~ 1 



-2(— + 3 - 2 J D)m 2 y p 2 + p 4 l Eob 2 ; 0, 0)5 (p 2 ; m 2 ,,, m 2 ^) 



+P 4 (4^V _ p 2\ Bo(p 2. 0> o)B (p 2 ; m 2 w /i w , m 2 w /£ w ), 

and F(p 2 ,m 2 v , D) and G(^,P 2 , -D) contain the generic two-loop contributions 

,40 



(8.6) 



F{p 2 ,m 2 w ,D) = m 2 w 

x Ti 3/4/ (mL) + 



2(-^ - 118 + 187L> - 100L> 2 + 18D 3 )m 2 w - 9(7 - 4D)(3 - £>)p 2 



40 

4(— - 96 + 61L> - UD 2 )m 2 w + 9p 2 



T 2 3'4'(m 2 y^ 

2 \ 1 1 o//i™2 



-20(l/£> - 1)(4 -8D + D 2 )m 2 w (m 2 w +p 2 )T V2 , w (m 2 w ) + 18(4™^ 
-p 2 )\2(l - D)m 2 w + (3 - 2D)p 2 }T 12y4 ,{m 2 w ,m 2 w ) + 9(4m^ -p 2 ) 



x 



4(1 - D)m w + 4(3 - 2D)m w p -p T 112Z >v(m w ,m w ,m w 



+72(1 - D)m^(m 2 w + 2p 2 )T 123 '4'5'(m 2 w ,m 



w, 



+8(1 - £>)m^ 2(m^ +p 4 ) + (8 - D)m w p T V7!wv {m 



17) 



G^ w ,p 2 ,m 2 w ,D) = —{D - Z)T iy y{m 2 w )+T 2yA ,{m 2 w /i w ) - 2\(^- - 1 



-2(-^ + 3-2D)m 2 v p 2 +p 4 



+ (3 - 2Z})p 2 ]Ti 23 /4/(m 2 y , m 2 ^/^) + 1 

^ii23'4'(^vv, w 2 ^, m 2 ^/^) 
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Figure 8.3: Higgs-dependent contributions to the Z self-energy 



The functions f{p 2 ,m 2 v ,D) and F(p 2 , m^,, D) are independent of the gauge parameter 
The occurrence of in g(iwiV 2 , m wi D) and G(^vk,P 2 , m w> D) indicates that the 
contribution from the "charged current" graphs is gauge dependent. This feature is ex- 
pected since also in one-loop order the graphs containing H^-bosons give gauge dependent 
contributions to the photon self-energy (see for example |TT| ). 

As explained above the integrals occurring in ( |8.4| ) can be solved analytically for 
general values of £w Using the expressions given in [|13|] the dependence on the variables 
p 2 , mly, D and can be studied. 

For all results contributing to the photon self-energy we also calculated the longitudinal 
part. It vanishes separately for each set of graphs specified above and without restriction 
on the values of the gauge parameters. We also studied the limit p 2 — > for all results 
given above and checked that they give zero for all values of the gauge parameters and 
the dimension D. 

Next we focus on the self-energy of the Z-boson. We first note that there are contri- 
butions from 27 graphs corresponding to the ones considered for the photon self-energy, 
i.e. from the diagrams shown in fig. ^T2] where the in- and outgoing photon is substituted 
by a Z. In addition to these, five graphs occur containing the Higgs-field, three of which 



are tadpole graphs. These diagrams are depicted in fig. |875 . 

We begin with the diagrams which are of the same form as in the case of the photon. 
The transverse parts of the "neutral current" graphs are treated in precisely the same 
way as described above. The results have the same form and contain exactly the same 
integrals as the corresponding ones for the photon self-energy. Since we are only interested 
in the structure of the results we do not list them here explicitly. As above, these graphs 
yield gauge invariant contributions. 

For the 15 "charged current" graphs we obtain: 

ESSV) = tt-^ e4( n"fj! 2 ^ 2 „ 4 lf(p 2 ,m 2 w ,m 2 z ,D) 



(4tt) 4 (1 - D) 2 12m 2 w m 2 z s 4 w 



-A 
' B 

+9g(€w,P 2 ,rn 2 v ,D)+F(p 2 ,m 2 v ,m 2 z ,D) 

+9(p 2 -m 2 z )(2m 2 w -m 2 z -p 2 )^G(ti W} p 2 J m 2 Wl D)^ . (8.9) 

The gauge dependent functions g(£vy? P 2 ■, m w > D) and G{^wi P 2 , m wi D) are the same as for 
the photon self-energy. They were given in flS.6Q and ( |8.8| ), respectively. The expressions 
for the gauge independent functions f(p 2 , m^y, m 2 z , D) and F{p 2 , rriyy, m 2 z , D) read 

f(p 2 ,m 2 w ,m 2 z ,D) = 2[(l/D-l)(4-6D + D 2 )(9m A w + (m 2 w -m 2 z ) 2 ) +9m 2 w (m 2 w 
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-2(2 - D)p 2 ) A (m 2 w )B (p 2 ; 0, 0) + (1 - D)(Am 2 w - m 2 z )(2m 



IV 



+m 2 z ) 



2m 2 w + (7 - D)p 2 Blip 2 ; 0, 0) - 9(2m 2 w - p 2 



4(1 



+4(3 - 2D)m 2 w p 2 - p 4 B (p 2 ; 0, 0)B (p 2 ; m\ 



Wi m w) ? 



D)m\ 
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For m| = the functions f(p 2 1 m 2 Vl m 2 Zl D) and F (p 2 , , m 2 z , D) coincide with the 
corresponding ones used for the photon self-energy, i.e. f(p 2 ,m 2 v ,D) and F (p 2 , m 2 ^ , D) . 

The characteristic feature of the result ( ^.9[ ) is that all gauge dependent contributions 
involving two-loop integrals are proportional to (p 2 — m 2 z ). This can be understood either 
by noting that p 2 = m 2 z is the mass shell condition or, in the framework of the intrinsic 
pinch technique ||10|| , by observing that this factor is necessary to cancel a lowest-order 
Z-propagator. 

In contrast to the photon case the light fermion contributions to the Z self-energy 
depend on the Higgs mass. Whereas the contributions from the three tadpole graphs are 
compensated by the usual choice of the renormalization constant, the first two graphs 
in fig. give a p 2 -dependent contribution to the Z self-energy. We obtain for the first 
graph 



^ZZ,(2) 
J T,H 



(P 2 



e 4 (2-D)(9- 12s^ + 8i 



w) 



(4tt) 



D) 2 A8p 2 c 4 v s 



4 
W 



+m 2 z T 23 , i ,(m 2 H ) 



2m| 



m 



H 



(3-D)(p 2 -m 2 H ) + (2-D)m 2 z 
m 2 z + (3 - 2D)p 2 T 123 >4>(m 2 z ,m 2 H ) 



T 13 / 4 /(m, 



+m 2 z 



m H — m 



2 \2 



2p 2 (m 2 H + (3 - 2D)m 2 z J + p 4 T 1123/A/ (m 2 z ,m 2 z ,m 2 H 
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which is obviously gauge invariant. The subscript H indicates the dependence on the Higgs 
mass which appears in the scalar integrals and as coefficients m 2 H and mjj, respectively. 
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(a) 



(b) 

Figure 8.4: The light fermion contributions to the W self-energy 
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The result for the second graph in fig. 3^3 is derived from (|8.12j) by using the appropriate 
coupling of the ti-quark. 

The result for the three tadpole graphs reads 



^ZZ,(2) 
J T, tad 



e 4 (2 - Dfm\ 



I8T1 



13'4' 



IV, 



+ 



11- 



[m 



w 



m|) 2 



+ 9 



in 



Tiv*{m 2 z )\, (8.13) 



showing that they are gauge invariant. This is also true for each of the graphs separately. 

In the study of the 7Z-mixing self-energy YT< z ^ 2 \p 2 ) no new features appear. The 
diagrams are of the same type as for the photon, i.e. we have contributions from the 27 
graphs shown in fig. |8.2| with one of the external photons substituted by a Z-boson. Since 
the results have precisely the same structure as in the case of the photon we do not list 
them in detail. 

The graphs relevant for the W self-energy are shown in fig. ^T4|. As can be seen 
in fig. |3.4| a there is no natural way to make a subdivision into "neutral" and "charged 



current" graphs. Like in the case of the Z self-energy we get contributions from diagrams 



depending on the Higgs mass. These are shown in fig. |8.4j b. 

We first consider the 41 graphs listed in fig. |8.4j a. The result in 't Hooft-Feynman 
gauge, i.e. £ 7 = £z = £^ = 1, is given in the appendix. Inspection of the integrals 
appearing in ( JA.lp shows that the only integral type not present in the results of the 
neutral gauge bosons is Ti 2 3'4'5'(m|, m%/). Like the integrals considered above it can be 
solved analytically. The result is given in |0|. 

The result for general values of the gauge parameters is very lengthy and will not be 
presented here. It depends on all three gauge parameters £ 7 ,£z and £vv- Therefore in 
this case also the neutral gauge bosons yield a gauge dependent contribution. In analogy 
to the Z self-energy all gauge dependent terms involving the basic two-loop integrals are 
proportional to (p 2 — rriyy). The integrals appearing in this result are the same as in the 



case of the 't Hooft-Feynman gauge (see ( |A.1| )). 

The results for the Higgs-dependent graphs in fig. |S.4|b are of the same form as those 



obtained for the Z self-energy in ( |8.12| ) and (|S.13 ), where in (|S.12 ) mz has to be substi- 
tuted by raw- 



9 Conclusion 

Radiative corrections are necessary for the comparison between theory and precision 
measurements. Whereas the methods for performing one-loop calculations in massive 
gauge theories are well established, so far no complete treatment of irreducible two-loop 
corrections has been achieved. In the SM these corrections are generally expected to be 
small, but a detailed analysis is necessary to determine their actual value. 

In this paper we presented a technique for reducing two-loop self-energies to standard 
scalar integrals in general massive gauge theories. It makes use of the tensor structure of 
the two-loop integrals, the symmetries of the scalar integrals and certain integral relations. 
The results are valid for all values of the invariant momentum p 2 , the particle masses, 
the space-time dimension D and the gauge parameters. Features like gauge invariance or 
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transversality are displayed directly at the algebraic level. We explicitly verified Slavnov- 
Taylor identities valid for the neutral gauge boson system of the SM by calculating several 
thousand Feynman diagrams. 

As an application we treated the light fermion contributions to the gauge boson self- 
energies in the SM. For the self-energies of the neutral gauge bosons the contributions 
associated with the exchange of 7 and Z in the fermion loop are gauge invariant. The 
VF-exchange yields gauge dependent contributions. At the pole, i.e. at p 2 = m 2 z for the 
Z self-energy, the gauge dependence of the generic two-loop contributions vanishes. The 
photon self-energy gives zero for p 2 = as required by the Ward identity. The Z self- 
energy acquires a Higgs mass dependent contribution even in the light fermion case. We 
showed that it is gauge invariant. 

The W self-energy receives gauge dependent contributions from all intermediate vector 
bosons. As in the case of the Z self-energy the generic two-loop contributions are gauge 
invariant at the pole, i.e. for p 2 = m^. The Higgs-dependent graphs form a gauge 
invariant subset. 

All standard scalar integrals needed for the treatment of the light fermion contribu- 
tions to the gauge boson self-energies can be solved analytically. We performed these 
calculations and obtained expressions in terms of polylogarithmic functions. The results 
will be presented in a related paper El. 
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A Result for the W self-energy 



The contribution from the 41 graphs listed in fig. |8.4] a to the two-loop self-energy of 
the IV-boson is given by 
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where the 't Hooft-Feynman gauge was used. 



B Integral relations 

We list here some relations needed for the calculations performed in this paper. The 
reduction formula for the integral V 2 345 reads 
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From this relation the formula for the integral Y^z^ can be obtained by taking the 



derivative with respect to m\. 
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For the integral Y 23A we get 
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The integral Y" 2 34 is expressible through and integrals of simpler structure via 
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The expression for the integral T 23 45/ 5 '(m 2 , M 2 , M 2 ) in terms of integrals with fewer 
propagators reads 
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-4M 2 (m 2 - p 2 )T 235 ,(M 2 , M 2 ) + (8 - 3D)(m 2 - M 2 )(M 2 - p 2 )T 23A (m 2 , M 2 , M 2 ) 



+ 



D(m 4 M 2 + M 2 p 4 ) - (4 - D)(m 2 (M 4 -p 4 ) -p 2 (m 4 - M 4 ) - 2M 6 ) 



-8m 2 M 2 p 2 T 2U5/ (m 2 , M , M ) - 4M'(m 2 - pTT V2U (m 2 : M , M : 



'>}• 



(B-4) 
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